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Abstract 

We study a model of heat conduction with stochastic diffusion of energy. We 
obtain a dual particle process which describes the evolution of all the correlation 
functions. An exact expression for the covariance of the energy exhibits long- 
range correlations in the presence of a current. We discuss the formal connection 
of this model with the simple symmetric exclusion process. 



1 Introduction 

Simple systems of particles on a lattice have received considerable attention in the last 
years, as they are a testing ground for exploring the properties of far from equilibrium 
states, for which at present no general theory is available. 

Amongst the most studied models are the so called exclusion processes in which 
particles diffuse amongst empty sites on a lattice. The problem has been attacked 
with techniques of statistical mechanics and probability, and many analytic results 
are available [T^l [23 E], in particular for the simple symmetric exclusion process 
(SEP) 

One is also interested in the transport of continuous quantities, especially energy. 
Kipnis, Marchioro and Presutti (KMP) [TTJ have introduced a model of energy trans- 
port in which the energy of neighboring sites is randomly redistributed. This model 



*Department of Mathematics and Computer Science, Eindhoven University, P.O. Box 513 - 5600 

MB Eindhoven, The Netherlands, c.giardina@tue.nl 

^CNRS-ESPCI, rue Vauquclin 10, 75231 Paris, France, jorge@pmmh.espci.fr 

■^Mathematisch Instituut Universiteit Leiden, Niels Bohrweg 1, 2333 CA Leiden, The Netherlands, 

redig@math.leidenuniv.nl 



1 



has also been thoroughly studied, although analytic solutions - in particular for the 
steady-state correlation functions - are harder to obtain than in the SEP. Both the 
KMP and the SEP have an interesting (and rather exceptional) feature: the evolution 
of the f^-point correlation functions can be exactly mapped onto a diffusion equation 
for K particles. This so-called 'duality property' [15] is a powerful tool, which in the 
study of the SEP yields the complete ergodic theory (see [To]). 

Despite the resemblance between energy and particle transport, studying both 
models is not a redundant exercise, as the SEP and the KMP models also show 
intriguing differences [I] , mainly in their large-deviation functions. Perhaps the most 
striking is the fact that in the implementation of an 'additivity principle' of the 
large deviation functions of two subsystems, the density at the interface has to be 
maximized in one case, and minimized on the other, in order to obtain the correct 
result [6]. 

In this paper we shall study a family of models of energy transport which, we 
shall argue, are the natural counterpart of the SEP models. In each site % of the 
lattice there is a free particle with momentum Xj. Between any two neighboring sites 
+ 1) and for any small time interval there is a random exchange of momentum 
that leaves \x\ + xj +1 } invariant. The same transport terms were already considered 
in models of wave propagation in random media [2T], [JJ. They later appeared as the 
high-energy limit of a chain with deterministic dynamics [9], and were also used as a 
stochastic perturbation (mimicking nonlinearities) of an oscillator chain [3], [8]. 

In the present paper we construct for our energy diffusion model a dual process 
that expresses the evolution of the i^-point correlation functions of the kinetic en- 
ergies in terms of a process of K interacting random walkers. We also give a closed 
expression for the stationary covariance (x?; x?) which confirms the presence of long- 
range correlations in the non-equilibrium stationary state (as already found before in 
the SEP PJ). 

The energy diffusion model is clearly very close from the physical point of view 
to the KMP model. But at the same time, as we have mentioned, it can be viewed 
as the continuation of the SEP family, in the following sense: on the one hand, the 
SEP can naturally be generalized [IB] to processes in which up to n particles are 
allowed per site (the usual SEP has n = 1), and hopping rates are proportional to 
the number of particles at the departure and the number of 'holes' at the arrival site. 
On the other hand, the energy diffusion model can be generalized to having in each 
site i of the lattice m free particles with momenta and the random 

exchange between x^ a and x i+ i^ leaves ^ 7 { x i,7 + x i+i,j} invariant. We shall see 
below that the energy transport model with m particles per site formally corresponds 
to the continuation of the SEP family to negative occupation number: m = —n/2. In 
other words, the energy diffusion model as defined above is, formally, 'the SEP with 
— 1/2 particles per site'. 

The rest of our paper is organized as follows. In section [2] we define our energy 
diffusion model. In section [3] we construct the dual process of interacting random 
walkers and derive the basic corollaries of duality, namely, existence and uniqueness 
of the stationary measure and expression of the stationary correlation functions in 
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terms of absorption probabilities. In sections 4 and 5 we derive the exact stationary 
temperature profile and energy-energy correlation function. In section 6 we show 
that our model satisfies local equilibrium. In section 7 we show a formal connection 
between our model and the SEP family, showing that the energy diffusion model can 
be viewed as a bosonic version of the SEP. In the last section we discuss possible 
asymmetric extensions of our model. 



2 The model 

The model is defined as a stochastic process on iV-dimensional vectors (x±, . . . , xn) G 
R N which have to be interpreted as momenta associated to lattice sites {1, . . . , N}. 
Additionally, lattice sites 1 and N are in contact with a heat reservoir at temperature 
T L resp. T R . 

The process is defined by its generator L (acting on the core of C°° functions / 
with compact support) 

N-l 

L = L 1 + L N + J2 L ^ (2- 1 ) 

i=i 

with 



d 2 f df 

L >^ T ^- x 'oi (2 - 2) 



of df 

L "^ TK d^- XN dr N (2 ' 3) 

w= (/) (2 ' 4) 

This corresponds, in the language of Fokker-Planck equation (or master equation) to 
the following evolution equation for the time- dependent probability density p(x,t): 

where L* is the adjoint (in L 2 (dx)) of L, i.e., more explicitly, 

N-l 



L* — L* + L* N + L iji+ i 



i=l 



d 2 f d 



d 2 f d 
L * Nf = TR dx^ + d^ {XNf) 
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Let us first motivate the choice of this generator (12. ip . The Li and part is 
the generator of the usual Ornstein-Uhlenbeck process which represents thermalising 
noise corresponding to heat baths at temperatures Tl, resp. Tr. 

To explain the other part, consider the operator 



d d x 2 



In polar coordinates x = rcos9, y = rsin9, this operator simply reads 

d9 2 

which means that in the process (x(t),y(t)) corresponding to (12.61) . r(t) = r(0) and 
8(t) performs a Brownian motion on the interval [0, 2ir]. More precisely, from the Ito 
formula, it follows that the generator (12. 6p corresponds to the stochastic differential 
equations (in Ito sense) 

dx{t) = -x(t)dt + V2y(t)dB(t) 

dy(t) = -y(t)dt-V2x(t)dB(t) (2.7) 

where B(t) is standard Brownian motion. In other words r 2 (t) = x 2 (t) + y 2 (t) = r 2 (0) 
with probability one and the angular variable 6(t) = arctan(y(t) / x(t)) is a martingale. 

The bulk part °f ^ ne P rocess { x i(t), • • • i x N{t)) corresponds then to the 

stochastic differential equations 

d Xl (t) = -x!(t)dt + \/2x 2 {t)dB 1>2 (t) 

dxiit) = -2xi(t)dt + V2x l+1 (t)dB hl+1 (t) - yftxi-itydBi-^t) ie[2,N-l) 
dx N (t) = -x N (t)dt-V2x N ^(t)dB N _ hN (t) (2.8) 

where B i i+ i(t) (i — 1, . . . , N — 1) are independent Brownian motions and B 0) i(t) = 
B N>N+ i(t) = 0. In this process, the total kinetic energy * s conserved, i.e., 

^ixfit) = with probability one as can be seen easily from Ito's formula. 

The full process, i.e., the process with generator (12.11) (including the boundary 
terms Li,Ln) corresponds to the system of stochastic differential equations 

dxi{t) = -2xi(t)dt + V2x 2 (t)dB lt2 (t) + \fW L dW(t) 

dxi(t) = -2xi(t)dt + V2x i+ i(t)dB iii+1 (t) - V^Xi^dB^it) i e [2, N - 1] 
dx N {t) = -2x N (t)dt-V2x N _ 1 (t)dB N _ 1<N ( y t) + ^lF R dW'{t) (2.9) 

where W(t),W'(t) are two independent Brownian motions, independent of all the 
other B Li+1 (t). 
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3 Duality 



The main tool which considerably simplifies the analysis of this model is duality. First 
note that the equations for the evolution of correlation functions of degree n for the 
x process are closed, i.e., the time derivative of the expectation of a polynomial of 
degree n in the variables x\, . . . ,x n does not involve expectations of polynomials of 
higher order [9]. 

We now show that the evolution in time of well-chosen polynomials reduces to a 
Markovian evolution of their indices. If we interpret the indices of the polynomials 
as discrete (indistinguishable) particle configurations (i.e., specifying for each site 
i € {0, . . . , N + 1} the number of particles), then the evolution of the indices turns 
out to become a jump process that conserves the total number of particles. This 
jump process is called the dual process. A similar situation arises in the case of the 
SEP (which is self-dual) [IS], in the case of an infinite system of independent random 
walkers (where the Poisson polynomials have the self-dual property) [10] . and in the 
KMP model p]. 

We index our polynomials by a vector £ = (6, • • • , 6v+i), 6 6 N and introduce 
the notation x = (xo, £i, . . . , xn, xn+i), where the first and last components are fixed 
by xq = VTl and xn+i = VTr- The polynomial D(£,x) is then defined by 



N 



26 



= T?Tj? +1 II (2&-1)!! (3 ' 1} 

Only even powers of Xi need to be considered, since other stationary correlations 
vanish due to the invariance of the generator under the transformation a;, — > — Xj. 
We interpret £ G f2 = N^" 1 " 2 as prescribing the number of particles in each lattice site 
i G {0,...,iV + l}. 

In order to introduce the generator of the dual process, we define, for £ G Q, 
i,j G {0, . . . , iV + 1}, the configuration 6'"' to be the configuration obtained from £ by 
removing one particle at i and adding one particle at j. On the dual variables £ G Q, 
we then define the following generator. 

Op(0 ;= 2tMt lfi )-<p{£)] + 26(26 + IMC 1 ' 2 ) - 

N-1 

+ Yl (26(26-i + l)^- 1 ) - y(0] + 26(26 + i + - ¥>(£)]) 

8=2 

+ 26r(26r-i + l)^^ -1 ) ~ + 2^^(£ W ) - <^(0] (3.2) 

where ip : Q — > K is an arbitrary function of the finite particle configurations. 

In words, we can describe the process generated by L as follows: a configuration 
£ = (6; • • • 1 6v+i) represents K particles (or walkers) on {0, 1, . . . , N + 1} with K = 
Yli^ 1 6- The walkers can only jump to neighboring sites and are stuck when arriving 
to sites or N + 1. The rate at which there is a jump of a walker depends on how 
many walkers there are at neighboring sites. If we have 6 walkers at site i, 6-i 
walkers at site % — 1 and 6+i walkers at site i + 1 (for % = 2, . . . N — 1) then each of 
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the walkers at site i jumps to site i — 1 at rate 2(26-i + 1) and to site i + 1 at rate 
2(2£j+i + 1)- At the boundaries: each of the 6 walkers at site 1 is absorbed at site 
at rate 2 and it jumps to site 2 at rate 2(26 + 1); each of the 6v walkers at site N 
is absorbed at site N + 1 at rate 2 and it jumps to site N — 1 at rate 2(26v-i + 1). 
Note that this process conserves the number of particles, i.e., 



For a single particle, i.e., £ = 5i, the dual process is then = 5x(t) where X(t) 
is a continuous time simple symmetric nearest neighbor random walk jumping at rate 
2 and absorbed upon hitting or iV + 1. For two particles, i.e., £ = Si + 5j the 
particles perform independent symmetric nearest neigbor random walks at rate 2, 
except when they are sitting at neighboring sites. In that case, i.e., if the two walkers 
are at neigboring places, one of them jumps to the place of the other one at rate 6 
(and other jumps are still at rate 2). 

Remark 1. This attractive interaction between the dual walkers has to be compared 
with the repulsive (hard-core) interaction between the walkers in the SEP. 

In order to formulate our duality result and its consequences, we denote by 
expectation in the dual process (i.e., the process with generator (13. 2p ) starting from 
6 We now can formulate the duality result. 

Theorem 1. Let £(£) denote the process with generator ( 13. 2p and x(t) the process 
(x (t) , x\ (t) , x 2 (t) , . . . , XN(t) , x/v+i (t) ) where {x\ (t) , . . . , xn (t) ) is the process with gen- 
erator ( 12.11) . and where xo(t) = \ZT£, XN+i{t) = VTr- Then we have 



N+l N+l 



l£(OI = E^) = E&(°) = l£(o)l 



1=0 4 = 



E s [D{Z,x{t))]=Es[D(Z(t),x)} 



(3.5) 



Proof. Start from ( 12.11) . For % 



1 , . . . , N — 1 we have 



L i)i+1 D{Z,x) = Ttr i 




x 




(2^-1)!! (26+1-1 

x 2 ^ x 2 ^ 1 



, 2 ?i+l 

i+1 



1)!! 



26+i(26 + l) 



(26-1)!! (26+1-1)!! 



i+1 




which implies 



L iii+1 D(£,x) 



(2{ i+1 (2{ i + l)[D(e +1 > i ,x)-D({,x)} 
+ 26(26+i + l)[D(C' i+1 ,x) - D{S,x)]) 
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Furthermore 



( 




Mi 1 } 



26[B« 1 '°,i)-OK,i)] 



and 



fc£{iV} v Sft 7 



2^jv-2 r 2 ^jv 

r R W-ll^Tijif-^ 2 ^(2^-l)!! 
2^D(^ + \ *)-£(£,*)] 



Therefore, we obtain 

LD(£,x) x) (3-6) 

where in the lhs of (13.61) the operator L is working on the x variable, and in the rhs 
the operator £ is working on the £ variable. 

This relation then lifts to the semigroups and the processes via a standard argu- 
ment, see e.g. the proof of theorem 1.1 p. 363 in [15]. (Notice that informally, by 
(ED E x D(^x(t)) = e ti D(£,-)(f) = e t£ D(-,x)(£) = %D(£(t), £)). □ 



Notice that in the dual process all particles are eventually absorbed at one of the 
boundaries 0, N + 1. Hence, the limiting configuration of the dual process starting 
from £ G Q will always be of the form k5o + IS^+i where k + l = J2i We say that a 
function / : Q — > M is harmonic for the process generated by £ if /(£) = E^[/(£(t))]. 
Since the configurations of the type kd~ + IS^+i are absorbing, the corresponding 
absorption probabilities 

cw(0 = Pe(C(oo) = + (3.7) 

are harmonic. 

The following propositions gives some consequences of duality. 

Proposition 1. If Tl = Tr = T ; then the unique stationary measure of the process 
x(t) = (xi(t), . . . , arjv(t)) zs i/ie Gaussian product measure with probability density 
function 

1 / A x 2 \ 

Pt{x) = w^W ew \U^) (3 ' 9) 

From every initial condition [x\, . . . , X/v) the process converges in distribution to pr- 
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Proof. We first show stationarity of p?- If T L = T R = T, for the Gaussian measure 
we have 



where |£| = Y^i=id- Therefore, using ( 13. 5ft 



E 2 [D(^,f(t))]pr(a;)da; = E € J D(£(t), x)p T (x)dx 



T® = J D{£,x)p T (x)dx (3.10) 



where in the third equality we used that in the dual process the number of particles 
is conserved. Next, to prove the convergence, remember that in the dual process 
eventually all particles are absorbed at or at iV + 1. Therefore, if Tl = Tr = T, 

]imEg(D(£,x(t)) = J" P 5 (£(oo) = k5 + l5 N+1 )T k T l = 

k,l:k+l=\$\ 

D(t,x)p T (x)dx (3.11) 

□ 

Proposition 2. For all ^ T R , there is a unique stationary measure Pt l ,t r which 
has finite moments of every order, and for this measure one has 

p TLtTR (D(^x))= TtT l R c kl (0 (3.13) 

k,l :fc+H£l 

where c^(C) ore the absorption probabilities defined in (I3.7p . 

Proof. If Pt l ,t r is invariant and has finite moment of all order, then, by (13.51) we have 

J D(^x)p TLiTR (dx) = Km J E x (D(£,x(t)))p; TL)TR (dx) 

(lim E e (T»(^(t),x)))/x Ti)Tfl (dx) 
^ T^T%P^{ OO )=k6 + l6 N+1 ) (3.14) 

k,l :k+l=\£\ 

where in the last equation we used that eventually all particle are absorbed at one of 
the absorbing states or JV+1. Conversely, if (I3.13P holds, then, since the absorption 
probabilities are harmonic for the dual process, by (13.51) we have 

J D(Z,x)p TL ,T R (dx) = J %[D^(t),x)]p TL>TR (dx) = J E#[Z?(£, x(t))]pT L ,T R (dx) 

which shows invariance of Pt l ,t r - D 
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4 Temperature profile. 

The following proposition shows the convergence to and linearity of the limiting tem- 
perature profile (linearity of the stationary temperature profile has been shown before 
ill). 

Proposition 3. 

a) Starting from x G R , the local temperature at site i satisfies 

Mtf®) = %{x 2 x{t) ) = X>(i,J>3 (4.2) 

3 

where X(t) is continuous-time simple symmetric random walk jumping at rate 
2, and absorbed atO, N + \, and where pt(i, j) denotes the transition probability 
of this random walk to go from i to j in time t. 

b ) The stationary temperature profile is given by 

T,:=^ Tr (4) = T L P i (X(oo) = 0)+T R f i (X(oo) = N+l) 



1_ ivTlJ +Tr {it+t) <■**> 

c) For the stationary local temperature autocorrelation function we have 

N 

E^, Ta (x?(0)x?(t)) -TiTj = ]TP,(X t = k,r> t^T^xl) ~ W,(r > t) 

k=i 

(4-4) 

where r is the absorption time of the random walk {X t : t > 0}. As a conse- 
quence, 

E„ T T (xj(Q)x 2 Jt)) — T{Tj 1 i, 



F j (r>t) N k | 



Proof. 



a) This follows from (I3.5P and the fact that the dual process starting from a single 
particle at % is given by the random walk X(t). 

b) This follows by taking the limit t — > oo in (14. 2 p and the well-known formula for 
the absorption probabilities of the random walk X(t). 

c) The first statement (I4.4j) follows from (j4.2p and (14.31) . The second statement 
(14.51) follows from the fact that the quasistationary distribution of the random 
walk X(t) is the uniform measure on {1, . . . , N}, i.e., for any function / : 
{1 V} -5: 



1 - 

lim E t (f(X t )\r>t) = -J2f(j) 



□ 
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5 Stationary two-point correlation function 

We now study the two point function ^i Tl>Tr (x^x^). In the dual process we then have 
two particles starting at sites i and j of which we denote the positions at time t 
by (X t ,Y t ). From (13 . 21) we see that these two particles behave as two independent 
random walkers (with rate 2) as long as they are not at neighboring sites. When they 
are at neighbor sites then on the next step they prefer to be on same site compared 
to being further separated. 

From (I3.13P we infer 

VT L ,T R {^)) = (1 + 2M (t 2 l %{{Xoo, Foe} = {0}) + Tl PyftXoo, Foe} = {N + 1}) 

+ T L T R ¥ ij ({X oo ,Y oo } = {0,N + l})^ (5.1) 

Unfortunately, now the absorption probabilities cannot be obtained by elementary 
probabilistic considerations since the two random walkers of the dual process are in- 
teracting (attracting each other when they are at nearest neighbor sites). We therefore 
proceed by directly solving the linear equations for y it j = ht l ,t r {x1x])- We restrict 
ourself to the case 1 < i < j < N since iiT L ,T R {xfx^) is obviously symmetric. We 
further denote yi = Ti = HT L ,T R {xf). 

The yij then satisfy the following system of linear equations: 

3 Vi-i,i - 2 yu + 3 y iti+ x = j = i 

Ui-i,i+i + Vi+i,i+i + Vi,i + Vi,i+2 - 8 y^i+i = j = % + 1 

yi-i,i + Vi+i,j + Vi,j-i + Vi,j+i - 4 y i:j = j > i + 2 (5.2) 

with boundary conditions given by 

Vi,N+i = UiUN+i = ViTr i = 0, . . . , N + 1 

Voj = VoVj = ThVj j = 0, . . . , JV + 1 (5.3) 

Plugging in the ansatz 

= A + Bi + Cj + Dij 
we find that these equations have the following solution: if j > i then 

rp2 I • Tr - T L f T R ~T l \ . T L (T R - T L ) . . (T R - T L ) 2 

y ^~ ±L + 1 N+l \ L+ N + 3 J + 3 JV + 1 + 13 (7V + l)(iV + 3) 

(5.4) 

and for the diagonal i = j we have 



yu = 3 [Tl — 



(Tr - T L f 
L (JV + l)(iV + 3) 



„ . T R — T L f Tr-T L 
+ 6 1 - ■ - T L 



N+l V ^+3 

2 (t r -t l ) 2 

+ 61 (N + l)(N + 3) ^ 
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The solution can also be written as a quadratic form in Tl and T R and then the 
coefficients of the quadratic form add up to 1 outside of diagonal and to 3 on the 
diagonal (as can be seen directly from (15. ip ). Explicitly 



Vi,j 



1 



iV + 3 



J 



1-1- 



N + 3 



N + l 



1 - 



Tl + 



»(2+j) r2 
(JV+l)(JV + 3) R 



iV+ 1 



i(2+j) 



(iV + l)(iV + 3) 



TlTr 



(5.6) 



and 



+ 
+ 



1 - 



iV + 3 
i(2 + i) 



N+l 



(iV+l)(JV + 3) 



T 2 



(iV+l)(JV + 3) (N + l)(N + 3) 
i \ f * i 



1-1- 



iV + 3 



1 



iV + 1 



z(2 + z) 



(iV + l)(iV + 3) 



T L T R (5.7) 



Comparing this with (15. ip we can read of the expressions for the absorption prob- 
abilities. E.g., for i < j 



— Y 



0) 



iV + 3 



1 



J 



N + l 



(5i 



Remark that this probability is larger than in the case of independent random walkers 



where the expression would be (1 
attractive interaction of the walker 
at the same site. This effect is how< 

Finally, we compute the covariance 



]v^rj)(l — at+t)- This means that due to the 
attractive interaction of the walkers, they have a (small) preference to be absorbed 
at the same site. This effect is however negligible (of order 1/N) as N — > oo. 



VT L ,T R {x 2 ix]) - /iT L ,T fl (a:- )/iT L ,T fl Gr- 



and we obtain 



2i(N + l-j) 2 

A-iR - J-L) 



(JV + 3)(JV + 1) 



(5.9) 



+ 



(2N 2 + 8N + 3)T| + 6T L T R -37% . 2{T R - T L ) ((3 + 2N)T L + 3T R ) 



(N + l)(N + 3) 



(N + l)(N + 3) 



2N (T R -T L f 
(N + l) 2 (N + 3) 



(5.10) 



Remark 2. Notice that the covariance cy manifests the presence of long range cor- 
relations which are believed to be typical for non- equilibrium steady states /7gj /. The 
fact that the covariance is positive is due to the attractive character of the interaction 
between the dual walkers. This has to be contrasted with the case of the SEP model 
where the covariance is negative and the interaction between walkers is repulsive. 
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We notice that in cy appears the Green's function of the simple (continuous time) 
random walk X t , 

2i(N + l-j) 
G{hj)= (N + l) 

i.e., the expected total time spent at i starting at j before being absorbed at or 
N + l. Moreover, for < «i < a 2 < 1, 

lim Nc Nai:Na2 = 2(T R - T L ) 2 oti{l - a 2 ) (5.12) 

N— *oo 

From (15.121) . in the spirit of |14j . one expects that in the limit N ^ oo the random 
distributions 

1 - 



N . . 

converge jointly to a Gaussian random field with covariance 

C(aci, a 2 ) = 2{T L - Tr) 2 «i(1 - a 2 ) + 2T 2 (a 1 )5(a 1 - a 2 ) 

where 

T(a) = lim T\ aN \ = (1 — a)T L + aT R (5.13) 

N — >oo 

with \_aN\ denoting the integer part of aN. 



6 Local equilibrium 

From the computation of the stationary two point correlation, we infer in particular 
that for a G (0, 1) 

lim ht l ,t r {xIn\) = 3T » = 7=r^ / ^i^dx 

N -*°° ^J2-kT(ql) J 

while for each i G N, i > 

Hm » TL ,T R (x 2 laNl xl aNi+i ) = T\a) = (^== / x 2 e^ 2T ^dx 

where T a is defined in (15.131) . 

This suggests that "around each macroscopic point a" , which we associate to the 
micro-point \_aN\ , there is a Gaussian "local equilibrium" distribution with variance 
T(a). More precisely we give the following definition cf. [IT] . 

Definition 1. Let Tl andT R be fixed and let a G (0, 1). We say that local equilibrium 
holds if for all n G N, for all ki, . . . , k n G N and for all ix, . . . , i n G {1, . . . , iV} fixed, 
we have: 



N^o D [Y1 k ^lcN\+i v X J HT LjTR (dx) = II / 
^ 00 "' \/=l / 1=1 J 



X 



2k, 



pT( a ){dx) = [T(a)] 



J \z=i / i=r x 1 ; 

where pxia) denotes the measure of a centered guassian variable with variance T(a). 
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This is equivalent with the requirement that 



lim T [aNi (fi TLtTR ) = G(T(a)) (6.2) 

iV — >oo 

where T|q,jvj denotes spatial shift, G(o) denotes the product measure on IR N with 
marginals that are normally distributed with mean zero and variance a 2 , and where 
the limit is in the sense that expectations of polynomials of type D(£,x) converge to 
the corresponding expectations in the Gaussian measure. 

The following lemma shows that factorization of the absorption probabilities is 
sufficient for local equilibrium. 

Lemma 1. Let Tl and Tr be fixed and let a G (0, 1). Suppose that for all n e N ; 
for all ki, . . . , k n e N, for all ii, . . . , i n 6 {1, . . . , N} fixed and for all K,L e N with 
K + L = YS =x h 



j™ P Er= a Mi^m (^(00) = K6o + L5 



\K 



N+l, 



lim F [aNl (X(oo) = 0) K P La7VJ (X(oo) = N+1) L (6.4) 

N—*oo 



then local equilibrium holds. 

Proof. Combination of (13. 13j) with (16. 4p gives 



T R (dx) 

J2 lim Phvj (X(oo) = 0) M P [aN] (X(oo) = N + l) K T^Tj< 

K,M:K+M=k 1 +...+k n 

lim (P La7VJ (X(oo) = 0)T L + P L ^j(X(oo) = iV + 1)T R \ 

J— too V / 



fel+...+fc„ 

(Pi^ArifXfoo) = 0m + Pi„vifXfcx)) = iV+l)T R ! 

□ 



lim T, fcl +r +fc " = T(a) fcl+ - +fc " (6.5) 



In order to see heuristically why the factorization property for the absorption prob- 
abilities holds, it suffices to see that if we start n dual random walkers Xi(t), . . . , X n (t) 
from initial positions \_aN\ I = 1, . . . , n then we can couple them with n inde- 
pendent random walkers X[(t), . . . ,X' n (t) started at the same initial positions such 
that for all e > and for alH 6 {1, . . . , n} 

X t (t) -X[{t) < eVt (6.6) 

with probability close to one for t large enough. Indeed, if dual walker I is absorbed 
at 0, then his happens at a time of the order N 2 , and then, for large N, (I6.6P tells 
that the corresponding independent walker is at distance less than eN from at 
that time. Therefore, the probability that this independent walker is absorbed at 
iV + 1 is at most eN/(N + 1) < e. So the probability (in the coupling) that dual 
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walker X\ and corresponding independent walker X[ are absorbed at different points 
is less than e. Therefore, if we have the coupling with property (16. 6p . we have for all 

e h ...,e n e{o,N + i} 

lim (Hx 1 (oo) = 6 U ..., X n (oo) = Ofi - P'(X((oo) = B x , . . . , X'Jpo) = 0,)) = 

where P refers to the probability measure on path space for dual walkers starting at 
positions [aN\ + ii, . . . , [aN\ + i n , and P' refers to the probability for independent 
walkers starting at positions [aN\ + ii, . . . , [aN\ + i n . 

To see that such a coupling exists, we observe that in the dual process there 
is only interaction of the walkers when they are at neigboring positions. In that 
case they jump with higher rate (than independent walkers) to the same position. 
The coupling then consists in letting the independent walkers and the dual walkers 
perform the same jumps, and having extra jumps for the dual walkers (which are 
not performed by the independent walkers) when they are at neigboring positions. 
The total time that dual walkers are at neigboring positions in the time interval 
[0,t] is less than t? +s with probability close to one. So the difference between the 
position of the independent walker and the dual walker is a sum of the order of tz +s 
independent increments taking values ±1, which is bounded by (t5+ 5 )2+ 5 < tz~ 5 
with high probability. 

A similar idea of coupling has been implemented in the context of the simple 
symmetric exclusion process, see [5]. For the full proof of the factorization property 
along these lines we however refer to future work. 



7 The energy-diffusion model and the SEP family: 
bosons versus fermions. 



In this section we shall see the formal relation between the SEP family and the 
energy-diffusion model. This can help as a guide to see the similarities in methods 
for treating both cases. 

We shall first consider a generalization of the energy-diffusion model to the case 
in which there are m momenta per site, and kinetic energy is exchanged between any 
two momenta in neighboring sites. The generator is again (I2.ip . with now 



L 1 = £ 



a=l 



Tr 



d 2 



d 



a=l 
j m 

m is 



T, 



'9xl a 



d 2 







R dx 2 



XN.a 



N,, 



dx N , 



d 



d 



(7.1) 

(7.2) 
(7.3) 
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The factor 1/m multiplying L iyi+ i is rather arbitrary, it sets the time-scale. We wish 
to show the connection of this family of models (labelled by m) with the (partial) 
exclusion process, with maximal occupancy n, and in which the jumping rate is 
proportional both to the number of particles on the departure configuration, and to 
the number of 'holes' (n minus the number of particles) in the arrival configuration. 
The evolution operator of this process can be written as the Hamiltonian H of the 
spin j = n/2 ferro magnet [IB] 

H = -L* SEP (7.4) 

with 

L *sep = ) E ( J ^ J m + J i4 + i + 2 «Vi - 2j 2 ) (7.5) 

i 

+ a (J{ - J{ - j) + 7 ( Jt + J°i ~ J) + S(J£ -J° L - j) + P{Jt + J° L - j) 

The factor 1/j is analogous to the factor 1/m in (17. 3p . The operators J+, Jf, J° act 
on the Hilbert space corresponding to < r < n particles per site <S>j|r)j as follows: 

J+|r), = (2j-r)\r)i 
Jf\r)i = r\r)i 

J»i = (r-j)\r)i (7.6) 

They satisfy the commutation relations of the SU (2) algebra: 

[JlJt] = ±Jt ; [J-, /+] = -2 J° (7.7) 

and they can be transformed to the conventional SU(2) matrices by a similarity 
transformation. 

Representations are labeled by the squared angular momentum operator: 

J 2 \jM>=j(j + 1)\jM> (7.8) 

with j = n/2, so that the ordinary SEP (with (0,1) occupation) corresponds to a 
representation of spin 1/2. 

Going back to the generalized model defined above, we can rewrite the generator 
of the energy diffusion process as the Hamiltonian H = —L* with 



+ 2(T 1 K?-K°-^)+2(T L K+-K° L -j) (7.9) 



where we have defined, in each site: 

— ^ ^ ^ X i,a 
a 

d 2 



K° 



2 ^ dxj 



2 



t. a 



iJ^IjS-^ + ^t^-} (7 - 10) 



a 
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These satisfy the SU(1,1) relations: 

[K°,Kt]=±K± ; [K~,K+] = 2K° 



(7.11) 



(note the sign difference with respect to SU{2)). Representations are labeled in a 
manner analogous to SU{2): 



Kf\kM> = [K°f - ~[K+ Kr + K; K+]\kM >= k(k - l)\kM > 



K°\kM > 



m\kM > 



(7.12) 



To identify which is the representation (i.e., the value of k), we check the value of Kf 
as applied to the constant (which is the zero eigenvalue eigenvector of L*)\ 



m 



m ,m 



Kr\l>=0 ■ Ji?|l>=*|l>=-|l> ; Jtf|l>=*(*-l)|l>=-(--l)|l> 

(7.13) 

Hence k = m/4, and in particular k = 1/4 for the process with one velocity per site. 

Consider the coherent-state representation of vectors and operators (cf. [HI E21 
El): 



ip(zi) = < Zi\ip > with \zi >= e z 



10 > 



< Zi\ip > with \zi >= e^* Ji |0 > 



(7.14) 



(where |0 > is the state anihilated by J i or K i ). The group operators act on such 
states as: 



and 



< Zi\K?\if> > 

< Zi\K~\i) > 
< Zi \K°\iP> 

< Zi\J+\i) > 

< Zi\Jr\^ > 

<Zi\J°M> 



z1^- + 2k Zi 
dzi 

d i / 



< Zi\*p > 



9 

Zi-z h k 

OZi 



< Zi\lp > 



(7.15) 



2 9 

-z 2 i — + 2jz i 



d_ 



dz, L 

< Zi\lp > 



< Zi\lp > 







Zi a 3 

OZi 



< Zi\lj) > 



(7.16) 



Now, writing the generator of the energy diffusion model and the SEP models 
in these representations, it is easy to check that models with k = m/4 of one class 
formally map into models with —j = —n/2 of the other class, albeit with somewhat 
different boundary terms. 
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Another way to see the relation between these two families of models is to note 
that we can rewrite the generator L of the energy transport model with even m in 
terms of m/2 bosons at a and m/2 bosons bi a in each site 

(all other commutators vanish) and writing 

K t = E a k 6 L ; K i = E bia(lia 

a a 

K °i = \ E { a l a ™ + 6 £a.) + r ( 7 - 18 ) 

a 

We can similarly write the generalized SEP with allowed occupancy n with n fermions 
a ia and n fermions b ia in each site 

(all other anti-commutators vanish) and writing 

a a 

J " = lj2( a l a ^ + b\ a k a )-n (7.20) 

a 

Hence, SEP and energy diffusion models are essentially fermionic and bosonic 
counterparts of the same models. This is perhaps not surprising, in view of the 
fact that Poisson processes (as the SEP) are related to Fermions, just as Gaussian 
processes are related to bosons [2] . 



8 Asymmetric generalizations 

In this last section we briefly indicate how to introduce asymmetry in our model 
without breaking the bulk energy conservation law. This is in the spirit of introducing 
a bias to move in a preferred direction in the simple symmetric exclusion process. 

Consider as the first example two neighboring sites with momenta x, y. Going 
back to Eq. (12.61) . we may add to the diffusion term a bias: 

and its obvious generalization to a chain then becomes the asymmetric analog of the 
bulk part of fl2~T|) 

E L JK=E(^^7-^^) +E ^ x ^{ x ^r Xi+1 id (8 - 2) 

Possible choices for the drift function E(x,y) are: 
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1. E(x, y) = E is constant. In that case, the product of Gaussian measures (13.91) is 
still invariant if we add thermalising noise at both ends of the chain at the same 
temperature T. Moreover, in the case of different temperatures, the equations 
for the stationary correlation functions of energies at different sites are still 
closed, i.e., the equations for the expectations of polynomials of degree K only 
involve polynomials of degree K. However, duality is lost just as in the case of 
the asymmetric exclusion process. 

2. E(x,y) = E(x 2 +y 2 )- This corresponds to the most general situation where 
the product of Gaussian measures is invariant (under the circumstances of the 
previous item). In the case of unequal temperatures, equations for stationary 
correlations are no longer closed. 

3. E(x, y) = Exy. This form of the drift has the advantage of being a product of 
two .fT-operators in 1 17. 10p . which is also the case if one goes from the symmetric 
to the asymmetric exclusion process (in the quantum spin chain language). 
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